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Abstract 

In these lectures the introduction to algebraic aspects of Bethe 
Ansatz is given. The applications to the seminal spin 1/2 XXX model 
is discussed in detail and the generalization to higher spin as well 
as XXZ and lattice Sine-Gordon model are indicated. The origin of 
quantum groups and their appearance in CFT models is explained. 
The text can be considered as a guide to the research papers in this 
field. 



Preface 



During last few years I have given several short lecture courses on the ap- 
plication of (algebraic) Bethe Ansatz to the integrable models of quantum 
field theory in 1 + 1 dimensional space-time, see |], 0, ||. Being essentially 
similar in spirit they differ in the choice of particular topics and the techni- 
cal details. This course, which was delivered at the Institute of Theoretical 
Physics, University of New York at Stony Brook, is self-contained, and do not 
copy It can be considered as a short and non-technical introduction 

into the subject of integrable models. Interested students must turn to the 
vast literature. A minimal list of references is appended. A lot of references 
to historic as well as research papers could be found therein. 

The job of taking notes of the course and producing the lecture notes 
was undertaken by Kostas Skenderis. Without his help this text would not 
appear. I am very grateful to him for his devoted work. 

I would like to acknowledge the hospitality which was provided to me by 
Professor C.N. Yang in the Institute for Theoretical Physics, University of 
New York at Stony Brook. 
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1 Lecture 1 



1.1 Introduction. 

There are many different historical directions which led to the subject of ex- 
actly solvable model: (1) the study of magnetic chains initiated by Bethe in 
the early thirties, (2) the work of Onsanger and Baxter in 2 dimensional clas- 
sical statistical mechanics, (3) scattering theory in many body problem with 
factorizable S-matrix (Berezin, Yang). Here we will follow one more path 
which again leads to exactly solvable models, namely the inverse scattering 
method. 

Solitons are discovered almost a hundred years ago. These are particle- 
like solutions of non-linear differential equations. The prototype equation is 
the KdV (Korteweg-de Vries) equation 

Ut = uu x + u xxx , (1) 

where u t and u x denote differentiation with respect to t and x, respectively. 
It has solution which describes a localized wave which moves without dissipa- 
tion. The KdV equation resurrected in the early sixties through the work of 
Kruskal and others. They are the ones who named the particle-like solution 
soliton. 

What is very interesting is the existence of many-soliton solutions which 
opens new possibilities for the particle spectrum. In contrast with the usual 
field theory where to each field corresponds just one particle, here we have 
a novel situation where one field might generate more than one particle. 
This possibility might be familiar nowdays after the emergence of the string 
theory but back in the early seventies the soliton mechanism of emergence 
of mass spectrum was the first example of circumvention of the perturbative 
paradigm. 

Using the inverse scattering method, many different theories in 1 + 1 di- 
mensions were studied and solved, see for example monograph [Q. Among 
them the Bose Gas model, whose equation of motion is a nonlinear Schrodinger 
(NS) equation 

.dV d 2 ^ IT|2t 

,_ = __ + J |»|^, (2) 

where *&(x, t) is a complex field and g is the coupling constant, and the Sine- 
Gordon model. The latter describes a relativistic real field <p(x, t) whose 
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equation of motion is given by 



□0+ — sin/30 = 0. (3) 

The constants m 2 and f3 are a mass parameter and the coupling constant of 
the model, respectively. The solution of this model revealed a very interesting 
spectrum with the mass of the particles given by the formula 

16m .717 ,,, 

m n = sin — , (4) 

7 16 V ' 

where 7 = 8/(3 2 . These studies finally led to the so-called algebraic Bethe 
ansatz (1978). 

A natural question is what a 1 + 1 dimensional system has to do with the 
real world which is 4 dimensional. First of all, string theory teach us that 
2D physics (worldsheet) can be ultimately related with higher dimensional 
world (target space). In fact, string theories are described by conformal field 
theories (CFT) which can be viewed as ultraviolet limit of exactly solvable 
models. (Alternatively, one might think the massive exactly solvable models 
as deformations of CFT's). 

A second example where 1 + 1 physics enters comes from 4D QCD. There 
are evidences that the high energy scattering is mainly 2D phenomenon (Li- 
patov, Verlinde- Verlinde) . 

The similarities between non-linear sigma model in 2D and Yang-Mills 
theories in 4D (asymptotic freedom, mass through dimensional transmutation 
etc.) gives one more motivation for studying 2D models. One hopes that 2D 
physics can teaches us how to tackle corresponding 4D problems. 



1.2 The XXX model. 

Consider the one-dimensional quantum periodic chain with N sites. We as- 
sociate with each site n a local spin variable 




(5) 



where the components of the vector a are the Pauli matrices. The spin 
variables act on the Hilbert space h = C 2 . The full Hilbert space is the 
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tensor product of all local Hilbert spaces 

n = h m . (6) 

Its dimension is dim(7Y) = 2 N . The spin variable act on the full Hilbert space 
H as follows 

s n = I®I---®s®---I, (7) 

where / is the unit operator in the local Hilbert space h. The periodicity of 
the lattice means that s n+ N = s n . The Hamiltonian of the system is given 
by 

H = J2[(sn,s n+l )--}, (8) 

n=l 4 

where (s n , s n+ i) is the inner product between s n and s n+ i defined in the space 
of two spins, namely C 2 <g> C 2 . Notice that the Hamiltonian is manifestly 
negative. 

Notice that we have done all possible regularizations. The lattice spacing 
A and the finite volume N provides natural ultraviolet and infrared cut-offs 
in the theory whereas the choice of the quantum spin as our variable renders 
the Hilbert space finite dimensional. After we solve the theory we will wish 
to take the thermodynamic limit N — > oo. This is a non-trivial limit and it is 
related with the theory of infinite tensor product. However, in some cases of 
physical interest things are simpler. For example, in the case of ferromagnet 
the configuration consists mostly of, say, spin up configurations so only a 
finite number of spins could be turned down. 

We shall use the permutation operator 

P=^(I®I + a®a). (9) 
If we choose as a basis in C 2 the vectors 

i+>=m. i->=m oo, 



then the permutation operator is given by the matrix 

(11) 
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in the basis | + + >, | H — >,| h >,| > • The Hamiltonian can be 

written in terms of the permutation operator as 

1 N 

H=-J2(P n>n+1 -l). (12) 

Z n=l 

We now introduce the Lax operators L n a . Their significance will become 
clear later on. They are operator which act on the local space h n ® V (now 
we use the notation h n for a local quantum space assigned to the site n), 
where V is an auxiliary space. In our example h = C 2 is the same as V but 
in general it does not have to coincide with it. The L nA is given by 



L n ,aW = H J n® h) +i{s n ® B a ) (13) 
A + isl 



(A-I)/ + iP n>0 (14) 




(15) 

where = s 1 ± is 2 and A is complex parameter. The index n refers to the 
site position whereas the index a is always an auxiliary space index. The role 
of the parameter A will be made clear later. In the last line we treat L n a as 
an matrix on V with coefficients as operators in h n . 

The Lax operator can be interpreted as a connection on the ID lattice. 
Let 4> n be vectors from 7i <g> V assigned to each lattice site. We say that n 
is parallel if 

4> n+1 = -L n (f) n , (16) 

(where the factor 1/A is a suitable rescaling), so that the Lax operator just 
realizes parallel transport. If we take the formal classical continuous limit 
A — » 0, h — > (of course after we reinstate the h dependence) and define the 
continuous spin variable as 

s(x) = lim x = nA, (17) 



then using flTJD and (|lq) we get 



= lim I^±i^ = (18) 
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where S(x) = (s(x),a). But this is the equation which was used in the 
inverse scattering method applied to magnet by Takhajan|Q. So ([16]) is a 
natural lattice and quantum deformation of (|TE|). 

The next task is to establish commutation relations between the Lax 
operators. First we introduce our notation 

A-ikB mn -- (A ® -B)j m |fc n = (v4 B )i mj fcn ; 

where A 1 = A® I and B 2 = I £g> B. Since L na is a 4 x 4 matrix there are 16 
matrix elements whose commutation relations we want to study. All these 
commutation relations can be compactly written as 



R-ai,a2\^ AO-^n.ai (A)-£/n,a2 (aO — L na2 (/i)Zy nai (A)i? aijC[2 (A-/i). 



(19) 



This is an equation in V\ <g> V 2 <S> h n . The indices a x and a 2 and the variables 
A and [i are associated with the auxiliary spaces V\ and V 2 , respectively. The 
matrix R ai ,a 2 governs the commutation relation of the matrix elements of 
Lax operators. In our case it is given by 



R, 



0.1,12 



(A) 



A + i 



{\I ai ,a 2 +iP ai 



(12 , 



(20) 



There is a nice graphical representation of (|19|). The Lax operator can be 
represented as a cross of two lines (fig. 1), one which represents the auxiliary 
space index (thin line) and one which represent the site index (thick line). 
Then the graphical representation of (19) is given in figure 2. The right hand 
of ( |T9] ) is obtained from the left hand side by just shifting the n-line to the 
other side of the cross point of the remaining lines. 




n 




Figure 1. Graphical representa- Figure 2. Graphical representation of 
tion of the Lax operator L n ai . Fundamental Commutation Relations 

(IT9D. 
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Having interpreted the Lax operator as a kind of connection along the 
chain, we can introduce an operator which describes parallel transport once 
around the chain. This operator is the monodromy matrix T a (X) 

T a (X) = L N>a (X)L N ^ a (X) ■ ■ -L 1>a (X). (21) 

The monodromy matrix satisfies commutation relations identical with the 
ones satisfied by the Lax operators, namely 

(A - /x)T ai (A)T a2 (/x) = T a2 (n)T ai (X)R aua2 (X - ft). (22) 

These commutation relations can be easily proven by repeatedly use of (|19|). 
The graphical representation of ( p2|) (fig. 3) gives an alternative (and easier) 
proof of ( p2"D (the "train" argument). 




Figure 3. The "train" argument. 

We introduce a family of operators 

F(X) = Tr (T a (A)), (23) 

where the trace is over the auxiliary space. The operator -F(A) is acts on 
the full Hilbert space 7i (global operator). From its definition it follows that 
.F(A) is a polynomial in A of degree N. It follows from Q2"2"D that -F(A) and 
F(yu) commute 

[F(A),F(^)] = 0, (24) 

so there are N independent quantities which mutually commute. As we will 
see in a moment the Hamiltonian is among them. So we have integrals of 
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motion in involution for a system with N degrees of freedom. In the classical 
case it corresponds to the Liouville's definition of integrability. 

Consider the Lax operator at the specific point A = i/2. From ( |HD we 

get 

£n,a(A = i/2) = iP n ^ a . 



Hence, 



F(X = i/2) = i"Tr a (P N>a P N „ lia ■ ■ ■ P, 



La; 



Using P n , a P m , a = P m ,„P n ,a, Pn,m = P m ,n and Tr a P a = / we get 

F(X = i/2) = i N P 1>2 P 2>3 ---P N>N ^. 
But the right hand side is just the shift operator 

U = PipPzp ■ ■ ■ Pn,n-\) 
which has the property that, 

UX n = X n+ iU, 



(25) 
(26) 

(27) 
(28) 
(29) 



for every local operator X n . It follows that U is just the exponential of the 
momentum operator II 

U = expiU. (30) 

We now want to express the Hamiltonian in terms of the operators F(X). 
Since 

d 



dX 



we get 



dF 
~dX 



N 



;N-1 



^ Tr a (P/v,a • ■ ■ Pj,a ■ ■ ■ P\,a)i 



(31) 



\=i/2 j=l 

where the hat denotes omission. Comparing ( p7D and (|T2]) we see that the 
Hamiltonian can be written as 



\=i/2 



N 
~2' 



(32) 



Notice that since the P(A) commute for different A there is no problem with 
definition of the logarithmic derivative. 
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Let the monodromy matrix be 



A(\) B(\) 



T ^-[c\x) Dm)- (33 » 

We want to diagonalize 

F(X) = A(X) + D(X). (34) 

It is easy to show that 

N 

^=n®i+>« ( 35 ) 

n=l 

is an eigenvector of F(X). Indeed, since 

4|+>=0, (36) 

we get that C(X)Q = 0, and consequently that the matrix T(X)Q is a trian- 
gular matrix. Furthermore, using (|l5|) , we get 

A{X)Q = (A+Vo, (37) 

D(X)Q = (X--) N n (38) 
2 



so that 



F(X)n = [(x+ % -) N + (x- l -f]n. (39) 

Now we are going to show that the B(X) can be used as spectrum raising 
operator. This will be done in the next lecture. 
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2 Lecture 2 



In the last lecture we have seen that the vector f2 which consist of spin up 
configuration in all lattice sites is an eigenvector of F(X). The vector Q plays 
a role similar with the vacuum state of the harmonic oscillator. To generate 
other states we act with B(X) on Q. In that sense the operator B(X) is a 
kind of raising operator. However, the new vectors will be eigenvector of 
F(X) only for specific values of A's. 
Consider the state 

$({A})=S(Ai)S(A 2 )...S(A,)n. (40) 

To check if this it is an eigenvector of A(X) and D(X) we need the commu- 
tation relation between A(X), D(X) and B(X). These are obtained from (pT|) 
and are given by 

[3(A), £(//)] = 0, (41) 
A(X)B(ji) = a{X - fi)B{fi)A{X) + /3{X - n)B{X)A{fi), (42) 
D(X)B(fj) = 7(A - fi)B(p)D(X) + 5(X - fi)B(X)D(p), (43) 

where a(A) = (A — i)/X and 7(A) = (A + i)/X. The specific values of (3 
and 5 are irrelevant for our discussion. From ([IT]) we see that the $ is 
symmetric function of A's. If the second term in the right hand side of ([12]) 
and ( ]4"3] ) was absent then the vector $({A}) would have been an eigenvector 
of A(X) and D(X) for all values of A's. However, the presence of these terms 
generates extra unwanted terms when we commute the A(X) (or D(X)) to 
the right. For a particular set of A's all extra terms vanish and $({A}) is, 
at the end, an eigenvector of F(X). This set is determined by the solution 
of a transcendental equation. For this particular set of A's the eigenvalue 
equation reads 

(44) 

Since the left hand side of is an analytic function of A, so has to be the 



right hand side. However, the right hand side has poles at A = A;. All the 
poles are removed if the A's obey the equations 

A±iZV = ±=±±i k =i,...,i. (45) 

Ah — % L 1 _li Ah — Am — 1 
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In fact, this is exactly the condition for the cancellation of the unwanted 
extra terms mentioned above (one can find a formal proof in ||. Thus, 
self-consistency of ([0]) automatically provides the set of A's for which it is 
valid. 

Consider the shift operator 

"-'"-"fwU' («) 



We apply U to Since the second term in right hand side of (44) vanishes 
at A = i/2, the spectrum of U is multiplicative and, hence, the spectrum of 
the momentum operator II is additive 

n$ = £ p(A m )$, (47) 

m=l 

where 

p(A m ) = -ln- (48) 

z A m - i/2 

For A real the momentum ranges over [0, 27r] . The spectrum of the Hamilto- 
nian is also additive 

i 

H<f>=Y, h{\ m )$, (49) 

m=l 

where 

The parameter A can be interpreted as the "rapidity" . (Recall that the energy 
and the momentum of a relativistic particle can be parametrized in terms of 
the rapidity 9 as p = msinh^, E = mcosh.8.) Eliminating the parameter A 
we get the dispersion relation 

h (p) = ^( C0S P~ 1 )- ( 51 ) 

The momentum of the system is quantized due to the fact that we have 
our system in the finite box. This becomes manifest if we rewrite (|45| ) as 

e ip(X k )N = -q S (\ k -\ m ), (52) 

m=l 
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where 

S(X) = £±j , (53) 

is the two-particle scattering amplitude. (Compare with the free particle case 
where the quantization equation is expipN = 1.) The fact that the /-particle 
scattering amplitude in the right hand side of Q512| ) is expressed in terms of 
two-particle one is a manifestation of the integrability of the model. 

Since the Hamiltonian is negative, the particle spectrum does not de- 
scribe physical particles. However, if we reverse the sign of the Hamiltonian 
(ferromagnetic case) then we get physical particles which are called magnons. 

Another observable of our system is the global spin. The spin operator 
is defined as 

N 

£ = E s n . (54) 

n=l 

The spin operator appears in the 1/A expansion of the monodromy matrix 



T(A) = X N \l + ^-^ + 0(1)1 . (55) 
A A 



The O (1/A 2 ) term is related with Yangian symmetries 0. Consider now the 
commutation relations ([H|) in the limit A — > oo, \i fixed. It follows that 

[^7 a + E,T a (//)] = 0, (56) 

which means that the monodromy matrix is invariant under combined rota- 
tions in the full quantum and auxiliary space. One have in particular the 
following relations 

[S,F(A)]=0, (57) 
which implies rotational invariance of the Hamiltonian H, 

\E 3 ,B(X)) = -B(X), (58) 

which means that the operator B(X) turns down one spin and 

[Z + ,B(X)]=A(X)-D(X). (59) 

We evidently have 

N 

s 3 fi = —n, (60) 
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so that (B3) leads to 

N 

£ 3 $ = (--/)$, (61) 

where / is the number of B-operators acting on Q (see ([TOD). There are many 
comments in order here. First of all notice the difference between the odd 
and even case. When N is even we get integer spin states, whereas when 
it is odd we get semi-integer spin states. In particular the ground state can 
only appear in the even case. So contrary to the usual intuition, the nature 
of spectrum is different for the N even and iV odd case even in the limit 
iV — > oo. Next we can prove using (|59|) and the Bethe ansatz equations ( [i5|) 
that all $ states are highest weight states, i.e. 

S+0 = 0. (62) 

To get the rest of the states we act with £_. Since £_ commutes with 
the Hamiltonian H all the descendants have the same energy. Furthermore, 
equation (H) and (H) imply that I < N/2. 

One can easily check that fl44|) is invariant under complex conjugation. 
This means that for every complex solution A of (|44]) , A is also a solution. 
In addition all A's are different (Pauli principle) 0. This can be inferred by 
examining again the self-consistency of (f45|) . If two A's were equal then (fl5|) 
would have a double pole. So we would need (I + 1) equations in order to 
remove all poles which is an over-determined system (we have only I un- 
knowns). In fact one can prove that in the large N limit the (1 + 1) equations 
are incompatible. So the A's ought to be different. 

We want to study the vacuum state. From ( pTj ) we infer that / = N/2. 
In addition, more detailed investigation shows that all A's are real. We take 



the logarithm in both sides in (f4~5|) and we use a branch of the logarithm in 
the form ] A ' A 

- In ( ) = 7r — arctan — . (63) 

i v A — ia> a 

Equation (|45|) becomes 

N/2 

N arctan 2Aj = arctan(Aj — A/) + nQj, (64) 
i=i 

where Qj are quantum numbers which parametrize the A's. They are integers 
or half-integers. One can prove that A is a monotonic function of Qj so that 
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Qj — > <5™ ax in the limit A — > oo. In fact, for finite A?" the range of Q is 

_^ + I<Q<^_i, (65) 
4 2- V ~4 2' K ' 

so there N/2 Q's. This means that Q takes all the values from — Q max to 
Qmax, so actually we have Qj = j. The ground state consists of a completely 
filled Dirac sea. 

We now wish to take the thermodynamic limit N —>■ oo. We introduce a 
new variable x = j/N. In the thermodynamic limit Xj — > A(x). Assuming 
regularity of the function A(x) one can replace the sum in (]64]) by an integral 
within an (9(1/A^ 2 ) error: 



,•1/4 

arctan2A(s) = 7ra; + / arctan \X(x) — X(y) ]. (66) 

J -1/4, 



This a non-linear equation but it can be reduced to a linear one as we will 
now show. Differentiating equation ( |66|) w.r.t. x we get 

2X'(x) , . . . f 1 



1 + 4P 7r + A ' (a:) /l + (A( a :)-A(y))^ y - (6?) 
We introduce the density of states, p(A), defined as 

P(A) = (68) 

This is a natural definition since dx = p(X)dX. Using Q5"8|), equation ( |57D 
becomes a linear integral equation for the density of states p(X): 

2 f 00 1 

TT4^ = ^ (A) + L 1 + (A- / ,)^ ( " ) ^ (69) 

This equation was obtained by L. Hulthen can be solved by means of Fourier 
transform. 

The ground state energy in the thermodynamic limit is given by 
E = $>(A;) -> N J h{X{x))dx 

/CO 
h(X)p(X)dX. (70) 
-co 
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It can be shown that it is equal to 

E = -N\n2 = - \n2 N . (71) 

This value is related with the residual entropy at T = 0. 

Let us now consider excitations. In the iV-even case the first excitation ap- 
pears when I = N/2 — 1 (triplet state) and it consists of two vacancies (holes) 
in the Fermi sea. Each of them carries spin 1/2. The parameters Q are 
modified 

Qj 3 ■ "U ./,) • 0(.l .12). (72) 

where 9 is the Heaviside function. The parameters Qj have two jumps at the 
position of the holes. Equation (BB) is modified accordingly 



arctan2A(a;) = + tt [6 (J - ji) + 9(j - j 2 )\ + / dyarctan [\(x) - X(y)]. 

N J -1/4 

(73) 

and, hence, 



2 7r r°° I 

— 2 = + ^ (A - A,) + i(A - A,)] + rTI) _ ? ^) dft (74) 

where Aj, i — 1,2, are the rapidities for the holes. Thus, the density of 
states is also modified 

p(A; Ai, A 2 ) = p vac (A) + l(a(A - Ai) + cr(A - A 2 )) (75) 

where p vac (A) is the vacuum density of states and cr(A) is determined by the 
equation 

roo \ 

ttct(A) + y ^ 1+ (7(/i)d// = -tt5(A). (76) 

The energy of the system is the sum of the vacuum energy plus the contri- 
bution from the holes 

E = E vac + S(h)+S(l 2 ), (77) 
where _ 

h(n)<r(fji - A)rf/i. (78) 

-oo 
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Similarly, the relative momentum k(X) is given by 

/oo 
p(ji)<r(p - X)dfi. (79) 
-oo 

It ranges from to ix. The relative energy S(X) and momentum k(X) can 
be interpreted as the energy and momentum of the physical excitation. The 
dispersion relation is very simple 

£{k) = | sin A;. (80) 

The 1-particle excitation will come from the iV-odd sector. 

The scattering matrix S ex for the triplet excitation is determined by the 
equation 

In S cx ( A) = J In S(n)a{n — X)dfi, (81) 
where S(fi) is given in (0). It is given by 

Scx(A) = (82) 

where 

r(l/2 + zA/2) _ T(s) , S _I + ,A 

f[A) ~ r(i + zA/2) "r(i/2 + s )' S ~2 + Y [ * 6) 

It is interesting to note that the function |/(A)| 2 appears also in the repre- 
sentations of SL(2, R) as a Harish-Chandra factor. 

The singlet solution can be described in a similar fashion, but one make 
use of the complex solution of BA equations (see [Q). 

Up to now we have only considered the antiferromagnetic case (Hamilto- 
nian negative definite). Let us briefly consider now the ferromagnetic case 
(Hamiltonian positive definite). The physical vacuum coincide with the Fock 
space vacuum state (compare with the antiferromagnetic case where the phys- 
ical vacuum is constructed by filling the Dirac sea). In the infinite volume 
limit N —> oo equation (ESI) does not yield quantization of momentum. The 
energy (|5lD becomes positive after the change of sign of the Hamiltonian and 
corresponds to physical particles. The operators S + and £_ cease to exist 
and only 

N 

Q = Z 3 --, (84) 
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makes sense. Thus the SU(2) symmetry breaks down. The spectrum has 
also bound states with the dispersion law 

£m(p) = 2Af X + ~ x (1 - cosp), (85) 

where M is an eigenvalue of charge Q. 

Most of the techniques described so far apply to other integrable models 
as well. In general to define a model we need a group G, a representation 
p, and a spectral parameter. In our case the spectral parameter was defined 
on the complex plane. However, it may also be defined on a strip (non- 
degenerate rational case), or on a torus (elliptic case). We must introduce 
deformation parameters 7 and k to define these two cases. Thus one can 
classify all the integrable models in terms of the quartet 

(G,/o,7,k) 

For example the NS model corresponds to general spin, whereas the Sine- 
Gordon to general spin plus one deformation. 
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3 Lecture 3 

In this lecture we will consider higher spin generalizations of the spin- 1/2 
XXX model. Instead of spin- 1/2 variables at each lattice site we now have 
spin-S 1 variables. The local Hilbert space is modified accordingly 

h n = C 2S+1 . (86) 

One can define a Lax operator L n (X) in the same way as earlier, 

L ^-{ X tf A-k)' (87) 

but now acts on C 2S+1 (g> C 2 . The commutation relations between the L nja 's 
are given by 

Rai,a,2\^ ^)L nai (X)L na2 (/i) = L na2 [^L nai (X)R aia2 (A-/i), (88) 

where R aim is given again by 

c 

1/2 




Figure 4. Graphical representation of (]88f) . The thin lines represent 
"spin-1/2" lines, whereas the thick ones are "spin-S" lines. 

The monodromy matrix T(A) is defined similarly with the spin-1/2 case. 
The vector Q = Ti&>u n is again an eigenvector of the operator F(X) = 
Tr a T(A), where u n is a highest weight spin state 

s+u) n = 0; slu n = Su n . (89) 

Using again the same ansatz for the general eigenvector of F(X) 

mX}) = B(X 1 )B(X 2 )---B(X l )n, (90) 
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we get the Bethe equations 



^ Afc + iS ^jy _ \ k — \ m + i fa — i \ (9i) 

Ah — i i j Ak — Am — % 

If {A} is a solution of (|9lD then $({A}) is an eigenvector of F(X). However, it 
is not clear at all how to get a local Hamiltonian. Since the auxiliary space is 
different from the local Hilbert space we can not use the same tricks we used 
in the spin-1/2 case where we had L n>a = P n>a at some specific value of A. In 
order to circumvent the problem we change the auxiliary space to V = C 2S+1 . 
So, from a S 1 ® 1/2 representation we go to a S (8) S representation. The new 

O C M 1 O C M 1 

Lax operator L n j is defined on C £g> C and is called the fundamental 
Lax operator ||]. This Lax operator is graphically represented by a cross of 
two S-lines (thick lines). In order to find its explicit form we have to solve 
simultaneously the equations which are graphically represented by the figures 
5 and 6. 




Figure 5. The S-S-l/2 equation. 




Figure 6. The S-S-S equation. 

Solution to these equation exists and it is unique. The fundamental mon- 
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odromy matrix is defined as 



T f - L nJ ■ ■ -L n! f, 



(92) 



and is represented as in figure 7. 



Figure 7. Graphical representation of the fundamental monodromy ma- 



We wish to prove that the trace of the fundamental monodromy matrix 
commutes with the trace of the auxiliary monodromy matrix 



The first one will yield the local conservation laws whereas the second is 
the one which the Bethe ansatz gives. Having a solution of the equations 
described in figures 5 and 6, the "train argument" (fig. 8) immediately 
shows that ( |93"D is satisfied. 




Figure 8. The "train argument" for equation fl9"3|) . 

Let us describe the algebraic structure underlying the Bethe ansatz. We 
propose the following mnemonic rule. Let us suppose that we have an asso- 



trix T, 



[Tr f T f ,Tr a T a ]=0. 



(93) 
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ciative algebra A, and an operator R on A <8> A which satisfies a universal 
Yang-Baxter equation 

-R12-R13-R23 = R23R13R12, (94) 

defined on A® A<g> A. The indices indicate to which of -4,'s the operator 
R belongs. For example R13 means that the operator R belongs to the first 
and third A. Let p(s, A) be a representation of the algebra A characterized 
by a discrete index s (the spin in our case) and a complex variable A (the 
spectral variable). A representative example is the loop algebra where a 
representation is given by X a>n = T a X n , where T a are the spin matrices 
(however it is not our case). We define 

{p{s x , Ai) ® p(s 2 , X 2 ))R = R S ^(X U A 2 ). (95) 

We assume that R SlS2 (Xi, A2) depends only on the difference (Ai — A2). Con- 
sider the case h n = C 2S+1 , V = C 2 , and apply @ to 

p(a u X)p(a 2 ,p)p(n,a), (96) 

where the a-indices refer to the auxiliary space and the n to the local Hilbert 
space, as usual. Equation (0) becomes 

-Rai,a 2 (A — p)R ai ,n(X — a)R a2 . n (p — d) = R a2 ,n (p ~ Cr)R ai ,n ( A — a)R ai ,a 2 ( A — p) ■ 

(97) 

If we identify -R a ,n(A) with L n a (A) then ( ]9"7| ) is just equation (P5|). The value 
A = is assumed to be the specific point where L n a coincides with the 
permutation operator P n ^ a . We permute the indices in ([J|) according to the 
permutation 

(;-)■ 

This yields 

-R12-R32-R3I = -^31-^32-^12- (98) 

Acting in ( |9"S| ) with 

p(n 1 ,X)p(n 2 ,p)p(a,a), (99) 

we get 

-Rni,n 2 (A — p)R a ,n 2 {v — H)Ra,m(o- — A) = R a , ni (c ~ A)-R aj „ 2 ((T — p)R ni .n 2 ( A — /i) , 

(100) 
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or with the identification of R a , n (^) with L„ ja (A), 

L ni ,a(cr — ^) L n . 2: a(cr — fx) R nun2 (\ — fx) = R nun2 (\ — fx)L n2:a (a — fx)L nua (a — \) . 

(fOf) 

This is the exact form of the equation given graphically in figure 5. 

It is now straightforward to find the Hamiltonian of a spin-S 1 system. It 
is given by 

H = J2H n , n+1 = J2fs(Sn ■ Sn+l), (102) 
n n 

where ^ 

H ni! n 2 = ^"^^^ ra i' n 2 ("^) ^ = o^ ?ril ' n2 ' (103) 

and f s (s n ■ s n+ i) is polynomial of degree 2s. For example for a spin-1 system 
we have 

f 1 (x)=x-x 2 . (104) 

We now generalize our results to the case of the arbitrary complex spin 
variable S (Verma moduli). The spin matrices can be expressed in terms of 
Bose creation and annihilation operators ip* and ip which satisfy 

[r,ip] = -i, (105) 

as (Holstein-Primakov) 

4 = r n (2s-r n ^ n ) 1/2 , (io6) 

S- = (2S-r n ^n) 1/2 ^n, (107) 

s 3 n = Wn-S. (108) 

If ip*ip n is an integer then we are dealing with a finite subspace since neces- 
sarily 

< 25. (109) 

When S is complex s + , s~ and s 3 form a representation of SL(2, C) and not 
of SU(2). In the limit S — > oo the spin matrices limit to 

Jim 4 = V25<, (HO) 
Urns; = v 7 ^^, (111) 
lim 4 = (112) 
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Let us find the S — > oo limit of the Lax pair given in fl37|), 

■ LnW _ i ( X + is 3 n is- \ , . 

~ S[ is+ A - isl ^ (113) 



= j ( l\/S iJn/V2S \ 

+ ^ -rjV2s -ix/s )■ 1 ) 

Consider now the continuum limit An — > x. Let A = 1/5. The creation and 
annihilation operators are of order v^A- To see this note that the discretized 
version of a delta function is 

*(*-y)-% (115) 



so that ( [L05|) becomes 



VVi i $mn (116) 



'A VA A ' 
which proves our assertion. Defining the continuum ^-field as 



4,(x) = %, (117) 




' + M .,.1, ^ . ("8) 



equation ( |114| ) becomes 

iL(x, A) 



where we have absorbed a factor of % in the parameter A. But this exactly 
the Lax operator for nonlinear Schodinger (NS) equation]^]. Hence the NS 
model belongs to same class of model with the XXX magnetic chains. It is 
recovered in the S — > oo limit of the latter. The momentum and the energy 
in the NS model are given by k and k 2 , respectively. Let us see if we can 
get these values from the XXX model. We have seen that the momentum in 
ferromagnetic case is given by 

P (A) = iln(^f). (119) 
In the limit S — > oo, p(X) goes as X/S. Furthermore, the energy limits to 

*=*T*~§-£- (120) 
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Identifying k with an appropriately rescaled A (and subtracting the constant 
from ( 120 )) we see that the S — > oo limit of XXX model correctly reproduces 
the energy and momentum spectrum. Let us stress that here we differ from 
others, who prefer to view the NS model as a limit of the spin-1/2 XXZ 
model. 

The spin-1/2 XXZ model is defined by the Hamiltonian 



H ~ ( S n S n+l + S n S n+l + J S n S n+l) 



;i2i) 



where J is some constant different from 1. We will treat the XXZ model as 
a g-analogue of the XXX model. 

The notion of q- deformation was introduced by Gauss in the last century. 
He defined the g-deformation of x as 



x n 



q — q 
q-q- 



The limit q — > 1 gives back the original value 



limx„ 



x. 



If we introduce a new variable 7, 



then x„ is given in terms of the new variable 7 by 



X n 



sinh 7X 
sin 7 



One can now deform the Lax operator of the XXX model as 



sinh[7(A + is 3 )} is sin 7 



is + shi7 sinh[7(A — is )] 
The deformed R matrix looks as follows 



R 



( sinh[7(A + i)\ 

sinh 7A i sin 7 

i sin 7 sinh 7A 

sinh[ 7 (A + i)\ j 



;i22) 



(123) 



(124) 



(125) 



(126) 



(127) 
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Note that we have changed the normalization w.r.t. (pO]). If we impose the 
fundamental commutation relations fllTf ), we discover that the spin matrices 
obey deformed algebraic relations^], 

[s 3 ,s ± ]=±s ± (128) 

sin(2 7 ^ 
S1117 

This algebraic relations define what we now call quantum group. Actually, a 
better name would be deformed Lie algebra. In the mathematical literature 
the notion of deformation of algebraic structure is well defined. One starts by 
making the structure constants CV,... t-dependent, where t is a deformation 
parameter. The algebraic structure determined by the deformed structure 
constants is, in general, nonequivalent to the original one since some of the 
properties of the corresponding algebra are changed after deformation. A 
well known example is quantum mechanics. It has been proven that quan- 
tum mechanics is stable deformation of classical mechanics (the deformation 
variable is the h). In this sense, one can use the word "quantum" instead of 
the word "deformation" . Drin'field was the first who used the term "quan- 
tum" in this sense and thus coined the term "quantum group" |TD . 
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4 Lecture 4 



In the last lecture we have seen that the XXZ^ model can be viewed as a 
g-deformation of the XXX5 model. We have constructed the Lax operator 
L(X) of the XXZ model as a g-deformation of Lax operator of the XXX model 
and we have seen that it obeys the fundamental commutation relations, 



R{\ - /i)L x (A)L 2 (/i) = L 2 ^)L\X)R(X - //), 



(130) 



only when the spin matrices obey certain relations. These considerations led 
to the discovery of the quantum groups. 

Now we shall elaborate on this in more detail. It will be convenient to 



perform a similarity transformation in all entries in ([13 

L(X) ^Q(X)L(X)Q- 1 (X) 



and 



where 



R(X - fi) -> Q 1 (X)Q 2 (fi)R(X - ^Q\X)- l Q 2 {^)-\ 

-iX/2 



Q(A) 



q 







q 





iX/2 



(131) 
(132) 
(133) 



In matrix form Q 1 (X)Q 2 (fi) is given by 
Q\X)Q 2 ^) 



I ? -i(A+/i)/2 









i(A- M )/2 



q 






; (A-M)/2 









i(A+ At )/2 J 



(134) 



and only the middle block of it acts non-trivially in ( |132|) . So the new R- 
matrix depends on the difference (A — yu). Explicitly we have now 



L(X, 7 ) 



sinh[7(A + is 3 )} is shi7e 



is^ sin 7e 



sinh[7(A - is 3 )} 



(135) 



and, 



R 



I sinh[7(A + i)} 



V 





sinh 7A 

i sin 7e~ A7 






i sin 7e 
sinh 7 A 




At 







sinh[7(A + i)] J 



(136) 
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We introduce a new variable 

x = exp 7A. 

In terms of this variable the Lax operator decomposes as follows 

1 



xL. 



-L. 



x 



where 



7 H/2 



(q-q x )s 



L_ 



q 

-H/2 



-H/2 



q 

-(q-q 



Q 





-H/2 



(137) 



;i38) 



(139) 



(140) 



and we introduce a notation q 11 ! 2 = exp 27s 3 . Similarly, the i?-matrix can be 
written as 

(141) 



R(X) = xR + - -R. 

x 



where 



R4 



( q \ 

1 (q-q' 1 ) 

1 

V q ) 



and, 



R- 



PRl l P. 



(142) 



(143) 



P is the permutation operator. We now substitute ( |138| ) and (|141|) into 
( |130| ). In each side of the resulting equation seven different powers of x and 
y (y = exp 7/1) appear, namely 

{x 2 , y 2 , x 2 y 2 , x 2 /y 2 , y 2 /x 2 , l/(x 2 y 2 ), 1}. 

This implies seven different equations which relate R± with L±. So, for 
example, equating the coefficients of "1" we get 



R+L^L/ 2 , + R—L\L 2 _ = L/ 2 ,L^R+ + L 2 _L\_R—. 



(144) 
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The other six equations are given by 

R HT) L\Ll = L 2 ± L l ± R ±{T) (145) 
R + L\L 2 _ = L 2 _L\R + (146) 
R_L\L\ = L\L\R_, (147) 

where the notation R±(=) means that we have to consider all four cases, 
namely plus or minus in R, plus sign in both L's, and plus or minus in 
R, minus sign in both L's. So, equation (|145|) is, in fact, four equations. 
However, only three of the seven equations are independent. We choose as 
our independent equations the following three 

RL\L\ = L 2 + L\R, (148) 
RL l _L 2 _ = L 2 _L X _R, (149) 
R + L\L 2 _ = L 2 _L 1 + R +J (150) 

where R can be either R + or LL. Let us illustrate how we can derive the 
rest of the equations starting from the above three. In particular we derive 
the following equation 

R_L\L% = L\L\R_, (151) 

starting from the equation 

R+L\L\ = L 2 + L\R + . (152) 

We multiply both sides with the permutation operator P. The left hand side 
gives 

PR + L\L 2 + = R-}PL\L\ 

= RZ l L\L\P, (153) 

where in the first step we used ( |143| ). Similarly, the right hand side yields 

PL 2 + L\R + = L\L\R-_ X P. (154) 

Equations (|153| ) and ( |154| ) imply ( |151| ). All but equation ([[44]) can be derived 
similarly. Equation (|144|) can be checked using the so-called Hecke property 

R + -R- = (q- -)P, (155) 
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which can be checked explicitly from (|142 ) and ( 143|) . 



The relations ( |145| )- ([147D can be taken as the defining relations for the 
g-deformation of the Lie algebra SL(2) expressed in terms of its genera- 
tors q H ,s + and s_, entering in ( |139|) , ( |140| ) |TT| . They are homogeneously 
quadratic. The matrices R± play the role of the corresponding structure 
constants. The Yang-Baxter (Y-B) relation 

R^R l3 R 23 = RfR^R 1 *, (156) 

and analogous relations involving R_ guarantee that higher order relations 
follow from the quadratic ones. (One can derive ( |156| ) from fl9~T|) for n = 



in an exactly the same way as we find (|145|) - (|147|) from ( |130|) ). The relation 



( |156|) guarantees that ( |148j )-( |15Ti|) are the only relations defining the deformed 
algebra SL(2) q . Indeed, following the two paths in the diagram 

L 1 L 3 L 2 ^ L 3 L 1 L 2 

L l L 2 L 3 = y = L 3 L 2 L 3 . (157) 

L 2 L 1 L 3 ^ L 2 L 3 L 1 

we get the cubic relation 

R 123 L 3 L 2 L 1 (R 123 y 1 = R 321 L 3 L 2 L\R 321 )-\ (158) 

where R 123 and i? 321 are the left hand side and the right hand side of the 
relation ( |156|) , respectively. So, due to this relation the relation (|158|) is 
empty. Very general theory-categorical considerations show that if the cubic 
relations are absent (so if the Y-B relation ( |156j ) holds) then no new relations 
will occur in any order of L. 
Consider now the operator 

R = PR+. (159) 
The Yang-Baxter relation takes the form 

R12R23R12 = R23R12R23, (160) 

and the Hecke relation is 

R 2 = I + (q--)R. (161) 
Q 
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This is a g-deformation for the transposition generators of the symmetric 
group a lk , realized through the permutation operator P, i.e., 



P12P23P12 — PtrPmP'. 



23 -T12-T 23, 



(162) 



and 



Indeed, we have 



P 2 - T 



R 



q=l 



(163) 



(164) 



so that ( p.62| ) and ( |163|) are the q — > 1 limit of Q160Q and ( |161| ), respectively. 

The Y-B relations ( |156| ) are characteristic of the Braid group of Artin. 
The representation (|135|) which we associate with the Lie algebra SL(2) q can 
be generalized to SL(N) q . The R- matrix will be then a triangular matrix 
N 2 x N 2 and the Hecke relation still holds. Other classical Lie algebras lead 
to their own i?-matrix and some generalization of the Hecke relation. 

The SL(N) g-deformed Lie algebra is defined in terms of the correspond- 
ing i?-matrix by the relations (|145|) - (|147|) , where matrices of generators L± 
are triangular and the diagonal elements are arranged in such a way that 



/+/" 

l ii L (N-i),(N-i) 



(165) 



This gives the definition of the g-deformed Lie algebra SL(N) q (or 
corresponding to SL(N). One can give a similar definition of the g-deformed 
Lie algebras corresponding to all classical series B,C and D, see [ jyj| . 

After this general discussion we return to the SL(2) q algebra. We redefine 
R + by multiplying it with g _1//2 , 



R4 



fq 1/2 



V 





-1/2 










9" 



-1/2 



1/2 











q 1/2 J 



(166) 



and consider its blocks as the spin 1/2 representation of the matrix L + . We 
see that the g-deformation of the Pauli matrices are given by 



H/2 



{ l/2 

q 





-1/2 



(167) 
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and, 




168) 



(169) 



Let us note that here the structure constants (i?-matrix) appear in the fun- 
damental representation rather than in the adjoint one. 

Important property of quantum group is a comultiplication operation, 
which corresponds to the addition of spins in the q — 1 limit. Let L' ± and 
L" ± be two independent (commuting) set of generators. Then their matrix 
products L' ± L± satisfy the same relation as each of the L's.. Indeed, we have 



R{L' + L'[)\L' + L'[Y = (L' + L+) (L' + L+) R. 



(170) 



One can combine the generators of quantum group into one matrix instead 
of two, introducing a matrix L, 



L,L~ X . 



It satisfies the equation 



L} R L?R-. = R,^L 2 R+L X . 



(171) 



(172) 



This equation can be proven by using manipulations similar to the ones we 
have used earlier. Let us see, for example, how the left hand side arises. 
We start from equation ( |146| ) and multiply from the left with L^Rz 1 and 
from the right with Rz l {L}_)~ l (L 2 _)~ 1 R_ . Then the right hand side of ( |146j ) 
yields the left hand side of ( |172| ). The relation ( |171| ) is a kind of Gauss 



decomposition for the matrix L, since the diagonal of L is divided by ( |165|) . 

The classical limit q — > 1 or equivalently 7 — > is more transparent in 
terms of L. The R + matrix becomes for small 7 



where 



R + = I + i 1 r + + 0{ 1 2 ), 











^ 





1 

2 


2 











1 
2 





V 








\) 



(173) 



(174) 
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Similarly, one can define a matrix r_ from the expansion of R-. Then the 
difference between r + and r_ is a Casimir operator, 



C = r + 



T— = (7 



( 1 













-f 


2 








2 


-f 





V o 








1 J 



(175) 



We expand also the matrix L, 

L = / + 7 / + C( 7 2 ). 

Let us keep the Planck constant H, so that q 
commutation relations imply that 



H 2 - iH 1 



h[c, i 2 



(176) 

e llh . The fundamental 
(177) 



It is evident, that we got the relations of Lie algebra SL{2) written in 
terms of the "structure constant" C. The usual form 



[r, l b ] = ihe ab T, 



(178) 



is obtained if one introduces l a as 



Z = £A«. (179) 

a 

In terms of the universal enveloping algebra, generated by l a or L±, we can say 
that it was the comultiplication which was deformed and the multiplication, 
which was left intact. 

The opposite deformation takes place for the dual object, q deformation 
Lie group, corresponding to a given Lie algebra. Let us deformed it also in 
terms of matrix elements - coordinates on the group manifold - T = \\tij\\. 
The following relation 

RT i T 2 = T 2 T i R ^ ( lg0 ) 

make coordinates non-commutative, introducing a new non-commutative 
multiplication. The comultiplication 

T \T" -> T = T'T" \ (181) 
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is the same as in the classical case q = 1. 

It is possible to combine the generators L and T to define a bigger algebra. 
For that the following relation between T and L suffices 



T 2 l} = R.L 1 R~ l T' 2 . 



(182) 



which corresponds to the left action of Lie algebra on its Lie group. The 
system of equations ( |172| ), ( |180| ) and ( |182| ) defines the cotangent bundle 
(T*G) q . In the semi-classical limit we have 



CT 2 



(183) 



The set of relations (|172|) , (|180|) and (|182|) is covariant with respect to 
the left shifts 

T^TS, L -> SLS-\ (184) 
but then S is necessarily quantized, namely it obeys 

RS^ 2 = S 2 S 1 R. (185) 

Note the matrix elements of S commute with the ones of L and T. 
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5 Lecture 5 



Let us describe in more detail the system of relations (|172|) , ( |180| ) and ( |182j ) 
which defines a "quantum top" . As we will see later, it is the quantum top 
which drives the zero mode of the WZNW model. The phase space of the 
classical top consists of coordinates g which are group elements and momenta 
uo which are elements of the corresponding Lie algebra. The Lagrangian of 
the system is given by 

L = Trgg^ui - ^Trcj 2 , (186) 
and the basic Poisson brackets are 

{g\g 2 } = o, (187) 

{uj\uj 2 } = [C,u 2 ], (188) 
{u\g 2 } = Cg 2 , (189) 

where C = a ® a and, as before 

g 1 = g®I (190) 
g 2 = T®g etc. (191) 

Using the equations of motion 

u = (192) 
g = tug (193) 

we get for the time development of the system 

g(t) = exp[u(p)t]g(0). (194) 

In the quantum case our system of commutation relations from the Lec- 
ture 4 is 

L X R-_ X L 2 R_ = J R; 1 L 2 J R + L 1 , (195) 
RT X L 2 = T 2 T X R, (196) 
T 2 L X = R + L 1 RZ 1 T 2 . (197) 
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In the limit 7 — > we identify g with T and lu with I, where I is the first 
term in the 7 expansion of L, 

L = I + 7 l + .... (198) 

Then the quantum commutation relations ( |195| ), ( |196| ) and ( |197| ) limit to 
CIT87D , ([1881) and (|189|), respectively. (C = r+ - r_, see eq. (|lToD ). The time 
evolution of the system is described by the equation 

T{n) = L n T(0), (199) 

where the discrete index n plays the role of time. The interpretation of T(n) 
as a time-evolved coordinate is supported by the fact that the pair (T(n), L) 
obeys the same quantum commutation relations as the pair (T(0),L). 

Let us now consider the Wess-Zumino-Novikov-Witten (WZNW) model. 
The degree's of freedom are group elements g(x,t) and the left and right 
currents = d^gg' 1 and jjf = d^g^g, respectively. The model is de- 
fined on a cylinder R 1 ® S 1 or the sphere S 2 depending whether we are in 
Minkowski or Euclidean picture. The left and right currents separately gen- 
erate a Kac-Moody algebra. These are two independent algebras since left 
currents commute right currents. From now on we restrict our attention to 
the left sector. Let 

l=3a+3i- (200) 

Then the fact that / generates a Kac-Moody algebra is expressed through 
the following Poisson bracket between the /'s 

{l\x), l 2 (y)} = i[l 2 (y), C}5(x -y) + ^'(x - y), (201) 

where 7 is the coupling constant and is related to the level k c \ of the Kac- 
Moody algebra by 

k c \ = -. (202) 

7 

Here we write k c \ because quantum mechanically the level might renormalize. 

The quantum lattice picture for this looks as follows. Consider a chain 
with an operator L n attached to each site. Let i?±(n-m) be 

I J if ^ TfL 

R±(n- m) = < R± i{n = m (203) 
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We impose the following set of commutation relations for L n 

£n-R_( n _ m _l)£m-R-(n-m) = ^ + ( n _ m )^i2+(n-m+l)^- (204) 

It is easy to see that the shift of arguments in two i?'s generates a derivative 
of a delta function in the classical limit. Indeed for H — > 0, A — > 

(m—n) 

— > I + r y5 mn r±, (205) 

and, 

L n -> I + Al(x). (206) 
(Recall that the Lax operator is a kind of connection on the ID lattice, so 

L n ~exp / l(x)dx, (207) 
J A 

where exp / means a path ordered integral.) 
Furthermore, 

^5 mn ^6(x-y), (208) 

and, 

-^(S n ,m+i - $nm) -^S'(x-y). (209) 

Combining the terms of order A 2 in (|204|) we get (l20lD with the 5' term 
present. Quantum corrections make the following connection of deformation 
parameter q and Kac-Moody level k: 

%ix 

^ = ex P^' ( 21 °) 

Equation (|204|) is, in fact, three equations 

L l n L 2 n R. = R^L 2 n Ll (m = n), (211) 
L l n L 2 n+1 = L 2 n+l R + L l n (m = n+l), (212) 

[LlL 2 m } = for|m-n|>2. (213) 

It is amazing that the our system still "remembers" the Virasoro algebra. 
Indeed, let An be the quantum structure generated by the Lax operators 
L n in a lattice with iV lattice sites. Then one can consider a new quantum 
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structure An 1 generated by products of some neighboring L's, where N 1 < N. 
Then An> ^ An- In other words if, for example, Li, L 2 , L 3 , L 4 , L 5 obey 
the commutation relations ( |204[ ) so does the set L[ = LiL 2 , L' 2 = L 3 and 
L 3 = L^L 5 . Hence, the density of points in the lattice is irrelevant. This is a 
kind of reparametrization invariance in the lattice. 

Consider now the local field g(x) for the WZNW model. In terms of the 
holonomy 

j rx 

u(x) =exp / l(x)dx, (214) 
Jo 

and the corresponding holonomy for the right current v(x) it can be expressed 
as 

g(x)=v(x)g(0)u(x), (215) 

where g(0) is a zero mode of the local field. We want to find the lattice 
analogue of this expression. We introduce two vertex operator u n and v n , 

u n = L n L n _i ■ • -L x , (216) 

and an analogue expression for v n made from right sector L's. Then, 

g n = v n g(0)u n . (217) 

The vertex operator un is actually equal to the monodromy matrix Mn = 
LnLn-i ■ ■ ■ L\. It satisfies the fundamental commutation relations 

M 1 RZ 1 M 2 R^ = R- 1 M 2 R + M 1 . (218) 

which coincides with that of the g-deformed algebra (|195|) . The derivation 
follows from the algebra fl21ip -( [2T3] ). The commutation relations ( |218| ) does 
not depend on the lattice spacing A and presumably stay intact in the proper 
continuum limit. Thus the quantum Lie group naturally enters the conformal 
field theory as a monodromy of the local current. 

It can be shown now, that the local field g n is commutative 

9n9l = gl9n, n^m, (219) 

and periodic 

g n+N = g n , (220) 
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if the pair (g(0),M) constitutes the top, i.e. the following commutation 
relations hold 

Rg 1 (0)g 2 (0) = g 2 (0)g\0)R, (221) 

and 

g 2 (0)M 1 = R + M 1 RZ l g 2 (0). (222) 

This means, that the g-deformed top dynamical system is a complete sys- 
tem of zero modes of WZNW model and the spectrum of this latter field- 
theoretical models patterned by the spectrum of the former finite-dimensional. 
Indeed, the Hilbert space of the WZNW model can be written as 

TAvznw = ® ( 22 3) 

i 

which is a sum of tensor squares of the irreducible representation TCj of the 
Kac-Moody algebra of spin j, j = 0, 1/2, . . . k/2. An analogue formula holds 
for the Hilbert space of the top, namely 

^to P = E^® V ?'' ( 224 ) 

3 

where Vj are finite dimensional representations of the g-deformed Lie algebra 
SL{2) q . Moreover the embedding of SL{2) q into Kac-Moody described above 
allows to state, that 

Hj = Vj ® H , (225) 

thus separating the contribution of the zero modes and oscillator degrees of 
freedom. This formula gives natural definition of fusion rules 

Hi®Hj = Vi ® Vj ® H Q . (226) 

Note that the factor 7io appears only once. These fusion rules are based 
on a new comultiplication law between the Kac-Moody and the quantum 
algebra 1)1^] . 

Let us illustrate the relations among the various algebras generated by a 
single Lie algebra g with a diagram: 

g Level 

/ \ 
CS ~ KM k g q Level 1 

/ \ / \ 

? KM Kq g qjk Level 2 

Figure 9. Relations between the algebras generated by the Lie algebra g. 
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At the level 1 we have either an affine infinite dimensional algebra algebra 
KMk or a finite dimensional but deformed algebra g q . Both of them are 
parametrized by one parameter, namely the level k for the Kac-Moody alge- 
bra and the deformation parameter q for the quantum algebra. As we have 
seen these two algebras are intimately related. Furthermore, the Kac-Moody 
algebra is related with a Chern-Simon theory (CS). At the level 2 we can 
have a local Lie algebra in 2 dimensions, an infinite deformed Lie algebra, 
or a further deformed finite dimensional quantum algebra. For the last two 
options we have natural candidates, namely the deformed Kac-Moody al- 
gebra KMk tq and the elliptic Sklyanin algebra g q ^, respectively. However, 
the relation between them has not been completely understood yet|l3[. For 



the first option there exist no natural candidate yet. Maybe a double-loop 
algebra is a possible candidate but this not clear at the moment. We stop 
our speculations here. 

Let us now briefly show how massive models fit into our framework. We 
consider the Sine-Gordon model (SG) for concreteness. The Lax operator is 
the one of the XXZ model: 

L(\ 7 ) = _L_ ( sm h[7(^ + ^ 3 )] is - sin 7 \ 
^ ' ' sin7 y is + sin7 sinh[7(A — is 3 )] J' ^ ' 

However, the realization of s ± , s 3 is different. We parametrize them in terms 
of real canonical fields ir n and n , 

[vr n , <p n ] = -il. (228) 

The explicit form of the realization is the following 

1 e" i7r/2 (l + K 2 e- 2i ^)e- in/2 , (229) 



2k sin 7 

s + = (s~)\ (230) 
s 3 = <p. (231) 

It follows that the fundamental commutation relations are satisfied. In the 
limit A — > and with k = mA one recovers the Lax operator for the Sine- 
Gordon model with mass m. Therefore, the SG model really belongs to the 
same class with the XXZ models which in turn is a deformation of the XXX 
model. Here we again differ from others who have studied the SG starting 
from the XYZ model of spin 1/2. 
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However, if we try to repeat the same analysis as we did for the XXX or 
the XXZ model we run into a problem. The lower left matrix element of L n 
has no zero eigenvalue. To get around this problem we consider the product 
of two Lax operators 

L n L n -i = ^ q y) J ' (232) 

Indeed, L n L n _i becomes an upper diagonal matrix when applied to a suitable 
chosen vacuum state Q. Furthermore, the state Q is an eigenvector of A and 
D, 

AQ = a(A)0, (233) 
DVL = d(\)tt (234) 

The Bethe equation reads 



-a(X)\N/2 /cosh(A + uj - ry/2) cosh(A — u — i^/2)\N/2 
~d{X)) ~ V C0S h(A + uj + « 7 /2) cosh(A -u- i^/2y 



tt sinh(A -fi + h) , 236 > 
J yf sinh(A — ji — 27) ' 



where uo satisfies the equation 



1 



2cosh2cu = m 2 + — . (237) 



m 

Let us compare this result with the Bethe equation for the XXZ model. There 
the left hand side has the form 

sinh(A + ijS/2) 
sinh(A - i^S/2) 



The hyperbolic cosines in ( |235| ) cause no problem because we can shift the A's 
to convert them to hyperbolic sines. The right hand side remains unchanged 
since it depends only on differences of A's. We conclude that the Sine-Gordon 
model can be viewed as inhomogeneous spin (-1/2) chain. 

Up to now we have consider only discretization in the space direction. 
One can also consider discretization in the time direction. There we have 
light-like Lax operators L which connect the lattice points in the light-like 
direction (Fig. 10). 
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L ^ u ) 7 (n+l) 

space direction 




L(A-w)\ / L{\-u) 

light-like directions 



Figure 10. Discretization in the time direction. 



For further discussion we refer to the article |L4 . 

This concludes these lectures. We hope, that we were able to show, 
that the magnetic chains and the Bethe Ansatz contain a lot of potential 
for classifying and solving the integrable models of quantum field theory. 
The parameters, entering the description of magnetic chain, i.e. group G, its 
representation and one (or two) anisotropy parameters can be identified with 
dynamical field variables, coupling constants and mass parameter. All known 
integrable models could be put into this scheme. The models of conformal 
field theory then appear as a particular (massless) limit. 
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